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We present measurements of the branching ratio and of the polarization amplitudes in charmless 
B s — > <p<p decays using data corresponding to 2.9 fb _1 of integrated luminosity, collected by the 
CDF experiment at the Tevatron. The branching ratio in B s — > <j)<j) decays is measured relative to 
the normalization mode B s -» J/y/<j> to be 88{B S -> <j><j>)/&(B s -> J/y/<j>) = [1.78±0.14(s*af) ± 
0.20(syst)] x 10~ 2 . Using the experimental value of $(B S —> J/yf<j>) we determine the B s — > tj)y) 
branching ratio 

@(B S -4 00) = 2.40 ± 0.21 (stat)± 0.27 (syst) ±0.82(BR)} x 10~ 5 . 

The polarization fractions are measured for the first time in this analysis and found to be: 

|A | 2 = 0.348 ± 0.041 (stat) ± 0.021 (syst) 
|A|| | 2 = 0.287 ± 0.043(sfof ) ±0.011 (syst) 



0.365 ±0.044(stat) ±0.027 (syst). 
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1. Introduction 

B s mesons were initially studied by the LEP experiments and then by the CLEO experiment at 
T(55). More recently the KEKB accelerator has been running at T(55) resonance as well, enabling 
the Belle experiment to do B s physics. The largest B s samples, however, are collected by the CDF 
and DO experiments at the Fermilab Tevatron. To date, the Tevatron has delivered about 8 fb _1 of 
integrated luminosity while each of the two Tevatron experiments has recorded close to 7 fb _1 of 
integrated luminosity on tape. The most recent Tevatron results in B s physics include studies of 
rare decays [|l|] like B s — > }X}X, B s — > e]X or B s — >• §)X)X, CP violation in B s — > //V0 decays ^ ||] 
and CP violation in inclusive semileptonic B decays [Q]. 

In this paper we focus on studies of charmless B s — >• <j><j> decays, performed by the CDF ex- 
periment at Fermilab. We present measurements of the branching ratio [||] and of the polarization 
fractions [^] in these decays using data corresponding to 2.9 fb _1 of integrated luminosity. 

Charmless B s decays are still to be fully understood. They offer the possibility to test our 
current theoretical understanding and represent promising ways to search for physics beyond the 
Standard Model (SM). The B s — > 00 decay is part of the so called B — > VV family in which the 
initial state B-meson is a pseudo-scalar (spin 0) and the final state VV contains two vector mesons 
(spin 1). In the particular decay of B s to 00, the final state is a CP eigenstate. Such decays can be 
used to measure the B s decay width difference (AT S ) and the phase responsible for CP violation in 
the interference between decays with and without mixing. To conserve the total angular momentum 
in B s — y 00 decays, the relative orbital angular momentum between the two mesons in the final 
state must be either 0, 1 or 2. In the angular momentum space, there are various bases which can be 
used to analyze decays of pseudo-scalars to two vector mesons, but any formalism involves three 
independent amplitudes for the three different polarizations of the decay products in the final state. 
Measuring the polarization fractions amounts to an important test of the corresponding theoretical 
predictions. 

Within the SM, the dominant process that contributes to the B s — > 00 decay is the b — > sss 
penguin digram shown in figure [l[ The same penguin amplitude appears in other B — > VV processes 
which exhibit significant discrepancies between the measured polarization fractions and the SM 
predictions. Explanations involving both new physics scenarios as well as newly accounted SM 
effects have been suggested to explain the observations. However, none of the existing scenarios is 
convincing enough. To solve this "polarization puzzle" it is important to study as many B — > VV 
decays as available. The first polarization analysis of B s — ^ 00 decays, performed by the CDF 
experiment is presented here together with an updated measurement of the B s — > 00 branching 
fraction. 

2. Measurement of the B s ; — >■ 00 Branching Ratio 

The B s — > 00 decay was first observed by the CDF experiment in 2005 [[7|] using a data sample 
corresponding to 180 pb _1 of integrated luminosity. The first measurement of the branching ratio 
2$(B S — > 00) was performed with 8 signal events and found to be [lA±0.6(stat.) ±0.6(syst.)] x 
10~ 5 . The analysis was updated in 2009 with a data sample corresponding to 2.9 fb _1 of integrated 
luminosity. The data were collected by a trigger which requires two tracks displaced with respect to 
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Figure 1: Main penguin diagram that contributes to the B s —^<j><j> decay. 

the primary vertex to enhance the contribution from long lived B mesons and suppress backgrounds. 
From the same dataset, B s — > J /yip decay are reconstructed as well and used as a normalization 
mode. This normalization mode was chosen because it has a topology similar to the B s — > decay 
and so, the measured branching ratio will be free of uncertainties from B s and B c i production cross 
sections, as it would not be using a similar Bj penguin decay (e. g. Bj — > §K*). 

The B s — > 00 decays are reconstructed from to two 0(1020) vector mesons where each 
meson is reconstructed from the decay — > K + K . Similarly, B s —tJ/yty decays are reconstructed 
from a J/y and a meson, where the J/y decays to /i + jU~ and the meson decays to K + K . 
Both B s decays described above lead to four particles in the final state and all four particles come 
from one potentially displaced vertex. The B s —> J /y<p is important on its own because it may 
improve the measurement of CP violation previously performed with a sample collected with a di- 
muon trigger The displaced track trigger may add about 25% more B s — > J /y<p events which 
are unique to this independent dataset. 

The events are selected according to an optimization procedure designed to maximize the ratio 
5/ y/(S + B), where S is the number of signal events and B is the number of background events 
under the mass signal peak. This figure of merit ensures minimal statistical uncertainty on the 
branching ratio measurement and it was verified to also optimize the uncertainty on the polarization 
fraction measurement described in the following section. The signal events are simulated while 
the background events are chosen from the B s mass sidebands. The variables used for the signal 
selection are chosen based on their discriminating power between signal and background. They 
are verified to be un-correlated and to exhibit good agreement between data and simulation. The 
most important variables used for the selection of both B s — > 00 and B s — > J /yip decays are the 
transverse decay length of the B vertex projected along the B transverse momentum, the impact 
parameter of the B meson, the quality of the four-track vertex fit and transverse momenta of final 
state particles. In particular, for the B s — > J/y<j) decays, one of the two muons from J/y is required 
to be identified by the CDF muon systems. 

Apart from the combinatorial background which is suppressed by the optimization proce- 
dure described above, other physics backgrounds are present in this analysis. These physics back- 
grounds come from real B decays which are misreconstructed as either B s — >■ 00 or B s — > J/ytp. 
In the case of B s — > J /y<p decays, the main background is B° — > J/yK*°, where K*° — > K + K~ . 
When the pion from K*° decay is identified as a kaon, the misreconstructed K*° falls in the mass 
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m(K*k'K*K~) [GeV/c s ] m(JA|> K*K1 [GeV/c 2 ] 

Figure 2: Left: 00 invariant mass. Right: J/y<p invariant mass. For both mass distributions physics 
background contributions are shown together with overlaid fits to the data. 



region. The background fraction fj/ vK *o = N(B° — s> J/\j/K*°)/N(B s — > 1S estimated using: 

f d 38{B° -» J/\jfK*°) 38(K*° -» A' + 7r~) £W°(y/yA0) 
fj/ V K*°- fs &( Bs ^j/yfi <%(<!> ^K+K~) e J M ( 

where £ J lv K *° (// y/0) is the trigger and selection efficiency of the B° — > J / \f/K*° decay recon- 
structed as B s -^J/Y<j> and e J /v<l> [ & the trigger and selection efficiency for e J ^^ , both determined 
using simulation, fj and f s are the production fractions of the Bj and B s mesons. The fraction 
fj/yK* i s found to be 0.0419 ± 0.0093. For the B s (j>(j) mode, the physics backgrounds come 
from B° -»• ((>K* Q -> K + K~K + n~ and B s -> K*°K*° -»■ K-% + K + %-. Using methods similar to 
equation 2.1 (see equations 2 and 3 in [^]) we find that the contribution of the B s — > k*°K*° is 
negligible and the contribution of the B° — >• 0,^*° — >■ K + K~K + n~ mode is about eight events. 

An important step in this analysis is to measure the signal yields of both B s — > J /iffy and 
B s ; — > 00. After applying the optimization procedure described above, the corresponding B s mass 
peaks are shown in figure ||. We find 17 66 ± 48 (st at.) B s — yj/yty signal events and 295 ±20(jta?.) 
fi s — > J/Y<p signal events. 

The branching ratio 38 of the decay B s —^00 normalized to the well known 38 {B s — > J/Y'P) 
can be evaluated using the following equation: 

39{B,^M) _ N H m/Y-^i) d$ T0T 

x TmH. . vv\ x Si x b mu 



^>(B s ^J/ ¥ <j>) N JM " &8{$^m " ef*r " m " 

where Nj/y^ are Nqq the numbers of fi v — > J/y(j) and fi s — > 00 signal events. £^{^ and 
are the combined trigger and selection efficiencies. The term £™ r accounts for the efficiency of 
identifying at least one of the muons in the muon detectors. Using the above ratio, the uncertainties 
in the production cross section of the B mesons cancel out and several systematic effects due to 
detector and trigger efficiencies cancel as well, allowing a reduced systematic uncertainty in the 
measurement of the branching ratio. 

The efficiencies for both B s . — > 00 and B s J/y(j) channels are obtained by taking the ratio 
between the number of simulated events that satisfy the trigger and selection criteria and the total 
number of generated events. The efficiency for the muon identification is determined in a differ- 
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ent way than the trigger and reconstruction efficiencies because the simulation does not account 
properly for muon acceptance and the corresponding uncertainties would not cancel in the ratio 
of efficiencies. The muon efficiency is determined as a function of the muon momentum and it is 
obtained by using inclusive J/y — > /i/i decays reconstructed in the same dataset where either one 
or both muons have been identified by the muon detectors. 

The most important systematic uncertainties in this analysis are listed here. The uncertainties 
in the number of signal events due to variations in the fit mass range that account for the possible 
presence of unidentified peaking background near the signal peak and uncertainty in the shape of 
the combinatorial background as well as the parameterization of the signal mass peak with a single 
Gaussian function instead of two Gaussians, uncertainties on the physics backgrounds coming from 
errors on the corresponding branching ratios, uncertainties on the muon efficiency, uncertainty 
on the ratio of the trigger and selection efficiencies due to poor knowledge of the polarization 
amplitudes and the decay width difference between the B s mass eigenstates. The total systematic 
uncertainty, excluding the error on the B s — > J /yty branching ratio is 11%. The final ratio of 
branching fractions is: 

_ [ 178±0 .14(tfaf.)±0.20(^jf.)] x 1(T 2 (2.3) 



Using the experimental value of the B s — > J/y<p branching ratio we obtain: 

3$(B S — > (j>(j>) = [2.40 ± 0.21 (stat.)± 0.27 (syst.)±0.82(BR)} x 10~ 5 (2.4) 

where the last uncertainty (BR) is the dominant contribution and comes from the error on the 
B s — > J/yQ branching ratio. We note that the world average for the 8§{B S — > J/yfQ) = (0.93 ± 
0.33) x 10~ 3 is based on a single CDF Run I measurement that assumed the ratio between the 
and B° d fragmentation fractions f s /fd = 0.40. The central value of SS(B S — > J/vQ) is scaled 
to reflect the current value of f s /f d = 0.110/0.399 = 0.28 [§]. Consequently, we use 88 (B s -> 
j/ ¥ (j)) = (1.35 ±0.46) x 10" 3 . 

This result is compatible with the initial observation with substantial improvement on 
the statistical uncertainty. The result is also compatible with recent theoretical calculations [^] 
and @]. 

3. Measurement of the Polarization Amplitudes in B s .— > 00 Decays 

As already pointed out in the Introduction, in the B s — >• <p<p decay, the dominant diagram is the 
b — > s penguin shown in figure ||. The same penguin amplitude is also relevant in in other processes 
which have shown deviations from the SM predictions. Such effects are the difference in the CP 
asymmetries in Bj — > K + n~ and B + — > K + n° and the potential difference between the sin(2j8) 
measurements in b — > sqq and b — > ccs B° decays. 

The decay amplitude in B s — > (fxp decays can be expressed in terms of three independent decay 
amplitudes, which correspond to the three possible relative angular momenta between the two <p 
vector mesons. In this analysis we use the helicity formalism in which the polarizations of the two 
vector mesons are either longitudinal with respect to the direction of motion Aq or transverse rela- 
tive to the direction of motion. There are two transverse amplitudes An and A± corresponding to the 
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two polarizations being parallel or perpendicular to each other. The fractions of these amplitudes 
can be measured from the analysis of the angular distributions of the final state particles (the decay 
products of the two mesons). 

Taking into account the V-A nature of the weak interaction and the helicity conservation in 
QCD, it is expected that the dominant amplitude is the longitudinal polarization while the transverse 



component is suppressed by a factor of my/mj [10]. This expectation is confirmed in tree-level 



dominated b ->• u transitions like B° -> p + p~ [|Tl|, @, B + ->■ p°,p + [|13|] and B + -> cop + 0, but 
it is not confirmed in B — > <f)K*, afe-s-i decay. In this decay, the transverse polarization fraction 



is about equal to the longitudinal polarization [15, 16, 17]. This unexpected result is known as 



the "polarization puzzle". Explanations involving either new physics [[T^, ^] or SM corrections 



including either penguin annihilation [(9J, |10J, gOJ or final state interactions Q21| , |22|, [23J, |25|] have 
been proposed. Recent work [|^] based on QCD predictions finds the longitudinal polarization 
fraction to be in excellent agreement with our measurement if the penguin annihilation amplitude 
is fitted to the B ->■ §K* data. 

We present the first measurement of the polarization fractions in B s — > 00 decays using 
2.9 fbr 1 of CDF data. As a cross check, we also perform the polarization fractions in B s — > J/\jf(j> 
decays. Both measurements are performed on the data samples selected for the branching ratio 
measurement described in section 2. For either of the two decays, we refer to the two vector 
mesons as V\ and V2 and to their decay products as final state particles Pi and P2 from V\ and P3 
and P4 from V%. 

We use the helicity formalism to describe B s — >■ 00 decays. The x' and x" axes are defined as 
the directions of the V\ and V% momenta in the rest frame of the B s meson, respectively. We define 
the angle 6\ (62) as the angle between the x' (x") axis and the Pi (P3) momentum vector, defined 
in the rest frame of V\ (V2). The <1> angle is defined as the angle between the decay planes of the 
two daughter particles. The three angles cb = (di,02, < i>) completely describe the directions of the 
final state particles. The distributions of these angles are used to separate the three amplitudes 
and determine their corresponding fractions. The probability distribution function (PDF) used 
to describe the helicity angular distribution for the signal is obtained by integrating out the time 
dependence. The differential decay rate as function of the helicity angles is given by: 

d 3 A(a>) 9 1 



deb 3271 W 



[F e (a>)+F (a>)} (3.1) 



where F e = ^[|A | 2 /i(«) + \A\\\ 2 f 2 ((o) + \A Q \\A\\\cos(8\\)f 5 {(b)], F = ^|A ± | 2 / 3 (S), W = 

^ + Here, Fl and Fh are the decays widths of the B s mass eigenstates, fj are functions 
of the helicity angles cb and 5n is a strong phase defined as 8u = arg(A* A\\). The decay widths Fi 
and F H are fixed to the world average. Although our trigger gives a non flat acceptance as a function 
of B s proper decay time, this time-integrated approach has been verified to give biases smaller that 
the statistical uncertainty of the polarization fraction measurement by using simulation and the 
measurement of the equivalent fractions in the J/y<p control sample. A similar formalism is used 
to describe B s —>J/ y0 decays in the transversity basis p^]. The data samples and the optimization 
procedures are the same as the ones in the branching ratio measurement described in section 2. 

The observables measured in this analysis are the polarization fractions |Ao| 2 and |An | 2 as well 
as the relative strong phase between them 81 1 . The measurement of these observables is performed 
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Figure 3: Detector angular acceptance projections for the helicity angles cos(0i), cos(02) and <t> used for 
the angular analysis of B s — >• <j>(j) decays. 

using an unbinned maximum likelihood fit using as event-by-event inputs the reconstructed mass 
of the B s candidate and the reconstructed helicity angles. The mass distribution is used in the fit 
to discriminate the signal from background. The angular distributions separate between the three 
polarization amplitudes. The signal mass distribution has a width of 20 MeV/c 2 for the B s — >■ 00 
and 10 MeV/c 2 for the B s — > J/\j/<j>- In both cases the signal is parameterized with two Gaussian 
functions with the same mean and different resolutions. The mass background distributions are 
described by exponential functions. The PDFs used to describe the helicity angular distributions for 
the signal are described in [Q]. The observed angular distributions in both helicity and trans versity 
bases are different than the expected theoretical distributions due to detector acceptance effects. 
The angular acceptance is determined using simulated signal events. The projections on the helicity 
angles 6i , 62 and <J> are shown in figure ||. The background angular distributions are determined 
from the B s mass sidebands. These distributions are parameterized with empirical functions. The 
<I> distribution is parameterized with a constant function and the angles d{ and d 2 are parameterized 
with functions of the form 1 +Bxcos 2 (0) where B is a parameter determined by the fit. Before 
performing the measurement of the polarization fractions in B s — > 00, several tests of the unbinned 
maximum likelihood fit are performed. The fit is tested on pseudo-experiments where no biases are 
found and the uncertainties are in the Gaussian regime. The polarization fractions are measured in 
B s — > J/y<p decays used as a control sample: 

|A | 2 = 0.534 ±0.019(sfaf.). |A||| = 0.220 ± 0. 025 (st at.). (3.2) 

In this case the polarization fractions are found to be in good agreement with previous CDF mea- 
surements from a di-muon sample [|7|]. Finally, samples of B s — > 00 are generated and passed 
through the full trigger and detector simulation and then through the analysis selection. The polar- 
izations are measured in these samples and good agreement with the generated values is found. 
Finally, we measure the polarization fractions in B s — >■ 00 decays: 

|A | 2 = 0.348 ± 0. 041 (st at.), |A||| = 0.287 ±0.043(rfof.)- (3.3) 

The measured strong phase is cos (Sit) = — 0.91^q 13- The fit projections onto the mass and helicity 
angles are shown in figure ^ which shows very good agreement between the data distributions and 
the fitting functions. 

The main systematic uncertainties on the B s — > 00 polarization fractions come from the depen- 
dence of the angular acceptance on the decay width difference AT S , uncertainties on the lifetimes 
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Figure 4: Fit projections for the mass component and the angular components in B s 



decays. 



of the heavy and light B s mass eigenstates T# and Tl and the potential KK s-wave contributions to 
the angular distributions. 

The final results, including systematic uncertainties are: 



|A | 2 = 0.348 ±0.041(af£tf.) ±0.021 (ayrf.) 
|A|| | 2 = 0.287 ± 0.043 (stat.) ±0.011 (syst.) 
\Aj_\ 2 = 0.365 ± 0.044(stat.) ± 0.027 (syst.) 
cos(8 l{ ) = -0mt° \l(stat.)±0.09(syst.) 
The longitudinal and transverse polarization fractions are: 

f L = 0.348 ± 0.041 (stat.)± 0.021 (syst.) 

fr =0.652 ±0.041 (stat.)± 0.021 (syst.) 



(3.4) 
(3.5) 
(3.6) 
(3.7) 

(3.8) 
(3.9) 



It is clear from this measurement that the SM expected amplitude hierarchy |An| 3> |A|| | ~ \A± \ 
is not valid in B s — > <p<p decays. Instead, the observed relation between the polaiization amplitudes 
is given by: |Aq| ~ |An| > |Aj_|, which is similar to the measurements for the b — > s penguin 



transition of B — > <pK* decays |[l_5j, g8|, |29J] which were the origin of the polarization puzzle. 
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We compare our results with various theoretical predictions of the polarization amplitudes. We 
find that our central values are consistent within the uncertainty ranges with the expectations of the 
QCD factorization while they are not in good agreement with the expectation of perturbative 
QCD [[K]] and QCD factorization [gj. 

4. Conclusions 

We have presented an updated measurement of the B s — > <p<p branching ratio using a data sam- 
ple corresponding to 2.9 fb _1 of integrated luminosity. Using the same data sample, we measured 
for the first time the polarization fractions in B s — > <p<p decays. The measured amplitudes confirm 
the previously observed polarization puzzle in certain B -» VV decays. 

Each of the two Tevatron experiments have currently accumulated about 7 fb _1 of data and 
expect 10 fb _1 by the end of the Tevatron running in 2011. With a sample three times as large, 
CDF will improve the statistical errors on the polarization amplitudes in B s — > by a factor of 
two and will attempt to measure the decay width difference Ar s in this mode. Further studies of 
rare B s decays and CP violation in the B s will be improved with more data. 
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